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1 Introduction 

The first half of this chapter describes the development in mathematical models of 
Brownian motion after Einstein's seminal papers [1] and current applications to opti- 
cal tweezers. This instrument of choice among single-molecule biophysicists is also 
an instrument of precision that requires an understanding of Brownian motion be- 
yond Einstein's. This is illustrated with some applications, current and potential, and 
it is shown how addition of a controlled forced motion on the nano-scale of the 
tweezed object's thermal motion can improve the calibration of the instrument in 
general, and make it possible also in complex surroundings. The second half of the 
present chapter, starting with Sect. 9, describes the co-evolution of biological motil- 
ity models with models of Brownian motion, including very recent results for how 
to derive cell-type-specific motility models from experimental cell trajectories. 



2 Einstein's Theory 

When Einstein in 1905 formulated the theory that quickly became known as his the- 
ory for Brownian motion, he did not know much about this motion' He was looking 
for observable consequences of what was then called the molecular-kinetic theory of 

^ Just how much he knew seems an open question that may never be answered [2]. 
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heat. So he was not concerned about the finer details of specific situations. In fact, 
apart from dated mathematical language, his papers on Brownian motion [1] remain 
paradigms for how to model the essence of a phenomenon with ease and transparency 
by leaving out everything that can possibly be left out. 
The simplest version of his theory, 

i(<) = (2i?)H(t) , (1) 

for the trajectory x{t) of a Brownian particle, here in one dimension and in the lan- 
guage of Langevin [3, 4], works so well also for real experimental situations that its 
extreme simplicity may be overlooked: No simplification of this theory is possible. 
The white noise ri{t) is the simplest possible: 

For all t, t', (Tj{t)) = and {rjitMt')) = S{t - t') . (2) 

When this noise is normalized as done here — as simple as possible — the dimensions 
of X and 77 require that a constant with dimension of diffusion coefficient appears 
where it does in (1). Equation (1) is mathematically equivalent to the diffusion equa- 
tion, introduced by Fick in 1857, in which the diffusion coefficient D is already 
defined, and that determines the factor 2D in (1). The new physics was in Ein- 
stein's assumption that Brownian particles also diffuse, and in his famous relation, 
the fluctuation-dissipation theorem 

D = A:bT/7o , (3) 

which relates their diffusion coefficient D and their Stokes' friction coefficient 70 via 
the Boltzmann energy k^T. It is derived by introducing a constant external force field 
in (1), and assuming Boltzmann statistics in equilibrium. For a spherical particle, 

70 = eirpiyR , (4) 

where p is the density of the fluid, v its kinematic viscosity, and R is the sphere's 
radius. 



3 The Einstein- Ornstein-Uhlenbeck Theory 

Details left out in the model described in (1-4) will be found missing, of course, if 
one looks in the right places. For example, the length of the trajectory x{t) is infinite 
for any finite time interval considered^. Ornstein and Uhlenbeck [5, 6] showed that 
this mathematical absurdity does not appear in Langevin's equation [3], 

Consider an interval of duration t. Split it into N intervals of duration At = t/N. In each of 
these, the mean squared displacement of the Brownian particle is 2D At. So on the average, 
the distance travelled in a time interval of duration At is proportional to (Aty^^ <x N~^^^. 
Consequently, the distance travelled in a time interval of duration t is proportional to t^^^ oc 
N^^"^ . Let A'' 00, and the infinite trajectory has been demonstrated. The proof can be 
made mathematically rigorous in the formalism of Wiener processes, e.g., which is just the 
mathematical theory of Brownian motion. 
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mx{t) = -7oi(t) + -Fthei-raal(0 , (5) 

where m is the inertial mass of the Brownian particle, and the force from the sur- 
rounding medium is written as a sum of two terms: Stokes friction, —70 x, and 
a random thermal force Ftheimai = {2kBTjoY^^ T]{t) with "white noise" statistical 
properties following from (2). The random motion resulting from (5) is known as the 
Omstein-Uhlenbeck process (OU-process). In the limit of vanishing m, Einstein's 
theory is recovered. Together, they make up the Einstein-Ornstein-Uhlenbeck theory 
of Brownian motion. 

The OU-process improves Einstein's simple model for Brownian motion by tak- 
ing the diffusing particle's inertial mass into account. As pointed out by Lorentz [7], 
however, this theory is physically correct only when the particle's density is much 
larger than the fluid's. When particle and fluid densities are comparable, as in the 
motion Brown observed, neither Einstein's theory nor the OU-process are consis- 
tent with hydrodynamics. This is seen from exact results by Stokes from 1851 and 
by Boussinesq from 1903 for the force on a sphere that moves with non-constant 
velocity, but vanishing Reynolds number, through an incompressible fluid. Hydro- 
dynamical effects that the OU-process ignores, are more important than the inertial 
effect of the particle's mass. These effects are the frequency-dependence of friction 
and the inertia of entrained fluid. Stokes obtained the friction coefficient, (4), for mo- 
tion with constant velocity [8]. Brownian motion is anything but that. Also, mass and 
momentum of the fluid entrained by a sphere doing rectilinear motion with constant 
velocity is infinite according to Stokes solution to Navier-Stokes equation [8, 9]. This 
gives a clue that entrained fluid matters, and the pattern of motion too. 

But since Einstein's theory explained experiments well, this hydrodynamical as- 
pect of Brownian motion did not demand attention. Not until computers made it 
possible to simulate molecular dynamics. 

4 Computer Simulations: More Realistic than Reality 

In 1964-66 Rahman simulated liquid Argon as a system of spheres that interacted 
with each other through a Lennard- Jones potential [10, 11]. He measured a number 
of properties of this simple liquid, including the velocity auto-correlation function 
4){t) = {v{t) ■ v{Q)), which showed an initial rapid decrease, followed by a slow ap- 
proach to zero from below, i.e., there was a negative long-time tail. Several attempts 
were made to explain his results theoretically, with mixed success. 

In the years 1967-1970 Alder and Wainwright simulated liquid Argon as a sys- 
tem of hard spheres and observed hydrodynamic patterns in the movement of spheres 
surrounding a given sphere, though all the spheres supposedly did Brownian mo- 
tion [12, 13]. Using a simple hydrodynamical dimension argument, and supporting 
its validity with numerical solutions to Navier-Stokes equations, they argued that 
the velocity auto-correlation function has a positive power-law tail, (f>{t) oc tr ^/"^ in 
three-dimensional space. This result is in conflict with the velocity auto-correlation 
function for the OU-process, which decreases exponentially, with characteristic time 
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m/7. But the 3/2 power-law tail agrees also with Alder and Wainright's simulation 
results for a simple liquid of hard spheres doing Brownian motion. 

This made theorists [14] remember Stokes' result from 1851 for the friction on a 
sphere that moves with non-constant velocity: There are actually two Stokes' laws, 
published in the same paper [8]. Einstein had used the simplest one, the one for 
movement with constant velocity, so the effect of accelerated motion is not accounted 
for in his theory. Nor is it in the Ornstein-Uhlenbeck theory. However, acceleration 
of a particle in a fluid also accelerates the fluid surrounding the particle, in a vortex 
ring (in three dimensions, and two vortices in two dimensions) that persists for long, 
disappearing only by broadening at a rate given by the kinematic viscosity [13]. 
In this way the fluid "remembers" past accelerations of the particle. This memory 
affects the friction on the particle at any given time in a manner that makes the 
dynamics of the particle depend on its past more than inertial mass can express. The 
result is an effective dynamical equation for the particle, Newton's Second Law with 
a memory kernel, as we shall see. 

5 Stokes Friction for a Sphere in Harmonic Rectilinear IMotion 

The friction coefficient that is relevant for a more correct description of Brownian 
motion, differs from the friction coefficient that most often is associated with Stokes' 
name, (4), but it is actually the main subject of reference [8]. Stokes was not ad- 
dressing the hydrodynamics of Brownian motion in 1851, but the hydrodynamics of 
an incompressible fluid surrounding a sphere that does rectilinear harmonic motion 
with no-slip boundary condition, at vanishing Reynolds number, and with the fluid 
at rest at infinity. The equations describing this motion are linear, however, and any 
trajectory of a particle can be written as a linear superposition of harmonic trajecto- 
ries, by virtue of Fourier analysis [15]. So the flow pattern around a sphere following 
any trajectory can be written as a superposition of flows around spheres in harmonic 
motion, as long as the condition of vanishing Reynolds number is satisfied by the 
arbitrary trajectory. It is for a Brownian particle's trajectory, so Stokes' result for 
harmonic motion is fundamental for the correct description of Brownian motion. 

In general, the instantaneous friction experienced by a rigid body that moves 
through a dense fluid like water, depends on the body's past motion, since the past 
motion determines the fluid's present motion. For a sphere performing rectilinear 
harmonic motion x{t; /) with cyclic frequency uj = 27r/ in an incompressible fluid 
and at vanishing Reynolds number. Stokes found the "frictional" force [8], [9, §24, 
Problem 5], 
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^fiiction(i; /) = -70 (^l + ^(^; /) 

3wpR^6{f) + ^7TpR'^^i{t;f) (6) 

= -7stokes(/)i(i;/) ; 

7Stokes(/) =70 + , (7) 

where only the term containing i{t; f) = —i2Trfx{t; /) dissipates energy, while the 
term containing x{t; f) = — (27r/)^x(i; /) is an inertial force from entrained fluid. 
The notation is the same as above: 70 is the friction coefficient of Stokes' law for 
rectilinear motion with constant velocity, (4). The penetration depth 5 characterizes 
the exponential decrease of the fluid's velocity field as function of distance from the 
oscillating sphere. It is frequency dependent, 

5{f) = (i^M)^ = R{U/f)-- , (8) 

and large compared to R for the frequencies we shall consider. For a sphere with 
diameter 2R =1.0 pm in water at room temperature where i/ — l.Q prn^/ps, = 
v/{itR^) = 1.3 MHz. 

Note that the mass of the entrained fluid, the coefficient to it in (6), becomes 
infinite in the limit of vanishing frequency /, i.e., the flow pattern around a sphere 
moving with constant velocity has infinite momentum, according to Stokes' steady- 
state solution to Navier-Stokes' equations. 



6 Beyond Einstein: Brownian Motion in a Fluid 

The friction on a sphere that, without rotating, follows an arbitrary trajectory x{t) 
with vanishing Reynolds number in an incompressible fluid that is at rest at infinity, 
is found by Fourier decomposition of x{t) to a superposition of rectilinear oscillatory 
motions Using (6) on these, gives 

Ffric,ion(/) = -7stokes(/)(-i27r/)i(/) , (9) 

which Fourier transforms back to [15], 

^friction(0 = -70 i (10) 

^^TipR^itit) - enpR^f^^ f dt'{t - t')-^/^i{t') . 
So the Langevin equation (5) is replaced by [16, 17] 



mx{t) = -Ffi-iction (0 + -Fextemal(0 + -Fthermal(i) , 



(11) 
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where Fextemai denotes all external forces on the sphere, such as gravity or optical 
tweezers, and -Flhermai denotes the random thermal force on the sphere from the sur- 
rounding fluid. 

Several authors have derived expressions for the thermal force using different 
arguments and finding the same result 

FthermaK/) = (2fcBr Re7s,okes (/) ) ^ Vif) (12) 

see overviews in [18, 19] ^. Briefly, Brownian motion in a fluid is the result of fluc- 
tuations in the fluid described by fluctuating hydrodynamics [9, Chapter XVII]^. In 
this theory one assumes that the random currents split up into systematic and random 
parts, the former obeying (Navier-)Stokes equation, the latter obeying a fluctuation- 
dissipation theorem. From this theory one derives the expression of the thermal force 
on a sphere in the fluid. 

Note that this description did not invoke a scenario of randomly moving molecules 
that bump into the micro-sphere and thus cause its Brownian motion. This scenario 
is correct for Brownian motion in a dilute gas. It is of great pedagogical value in un- 
dergraduate teaching. But it does not apply to fluids! The scientific literature shows 
that some undergraduates proceed to become scientists without realizing this limita- 
tion on the scenario's validity. However, the coarse-grained description that replaces 
a molecular description with a hydrodynamical one, is a very good approximation on 
the length- and time-scales of the thermal fluctuations that drive the Brownian mo- 
tion of a micron-sized sphere in a fluid. This is why fluctuating hydrodynamics [9, 
Chapter XVII] is formulated by a "stochastization" [20, Sect. 4.6] of Navier-Stokes 
equation, and not by coarse-graining Langevin equations for individual molecules in 
the fluid. The correct physical scenario to bear in mind is one of molecules squeezed 
together "shoulder-to-shoulder" in a manner that allows only collective motion, sim- 
ilar to that observed in a tightly packed crowd of people. 

Equations (10-12) constitute the accepted hydrodynamically correct theory for 
classical Brownian motion, i.e., in an incompressible fluid. It differs from Einstein's 
theory in a manner that matters in practise with the precision that optical tweezers 
have achieved recently [21, 22]; see Fig. 1. 

Power spectra of micro-spheres in optical traps can be measured with stochastic 
errors below 1% [21]. So the differences in Fig. 1 between Einstein's simple theory 
and the hydrodynamically correct theory for Brownian motion in a fluid can be ex- 
posed experimentally [21, 23, 24]. The form of the thermal force in (12), on the other 
hand, remains a theoretical result. It is not a controversial result, it is not questioned. 
But because it is a small effect, it has not yet been demonstrated experimentally. 

' Here we have written the frequency-dependent noise amplitude explicitly, and to this end 
introduced r;(/), the Fourier transform of a white noise r;(t), normalized as in (2). 
Readers familiar with the Green-Kubo theory of linear response to perturbations may ap- 
preciate fluctuating hydrodynamics as a case where the order of linearization and ''stochas- 
tization^ [20, Sect. 4.6] is a non-issue by virtue of the Reynolds number for thermal fluctu- 
ations. 
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Fig. 1. Thie power spectrum of Brownian motion in an optical trap according to Einstein's 
theory, PLoicntz , divided by the hydrodynamically correct power spectrum for the same mo- 
tion, Pnydro; see [21] for explicit expressions for the two spectra. Fully drawn line: Trap 
with Hooke's constant 3.8-10^^ pN/nm for a micro-sphere with diameter 1 fj,m. Dashed line: 
Hooke's constant 1.9- 10"^ pN/nm for a micro-sphere with diameter 0.5 ^m. At low time reso- 
lution, i.e., low frequency, the error vanishes. Einstein made an excellent approximation when 
he chose Stokes' law for constant velocity to characterize motion along a. fractal trajectory. 



7 Power-Law Tails 

In the absence of external forces, the position power spectrum of Brownian motion 
following from (10-12) is 

/l~2|\ 2fcBrRe7stokes(/) 
P /) oc {\X CX I I TJTT^ ■ (13) 

\m{2nfy + t27r/7stokes(/)r 

Here, the frequency-dependent numerator is the power spectrum of the thermal force 
in (11), while the denominator is given by the other terms in (11). The frequency- 
dependent friction coefficient, 7stokes(/), appears both in numerator and denomina- 
tor, and both appearances contribute, with opposite signs, to the t^^^^ power-law tail 
in the velocity auto-correlation function. 

By Wiener-Khintchine's theorem, the velocity auto-correlation function is 

(/)(t) = (i(t)x(O)) CX r d/e-'2-*/(2^/)2p(/) . (14) 
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At asymptotically large values of t, (f>{t) is given by -P(/)'s behavior at small values 
off, 

(2^/)2p(/) = 2D{1 - if/Uy/'- + OiflU)) . (15) 

Hence 

•^C*) = — + Oit-^^'') for t^oo , (16) 

quite different from the exponential decrease following from Einstein's simple the- 
ory, but not conceptually different from it [14, 16, 17]. 

Experimental evidence for this power-law tail remained sparse for years. Dy- 
namic light scattering offered promise of its observation, but only Boon and Boullier 
[25, 26] reported an experimental result of the magnitude predicted theoretically, 
with statistical errors about half the size of the signal. Paul and Pursey used pho- 
ton correlation dynamic laser light scattering to measure the time dependence of 
the mean squared displacement of polystyrene spheres with radius i? ~ 1.7 /im [27]. 
They found clear evidence for the expected f^'^/^-behavior (t^/^ in the mean squared 
displacement), but with an amplitude of only 74 ± 3% of that predicted theoreti- 
cally. They never found the reason why 26% of the theoretically expected amplitude 
was missing [28]. Ohbayashi, Kohno, and Utiyama [29] also used photon correlation 
spectroscopy, on a suspension of polystyrene spheres with radius 0.80 /im, and found 
agreement between the theoretical amplitude of the t~^^^ tail and their experimen- 
tal results which has 9-10% error bars. Their results also agree with the predicted 
significant temperature dependence. This convincing experiment thus supports the 
vaUdity of the theory (11). This is the current experimental status of the power-law 
tail of the velocity auto-correlation function of classical Brownian motion. 

Or was, when this chapter was written. But before it went into print, Ref. [24] ap- 
peared. Strangely, the velocity auto-correlation function is not given in [24], though 
its authors have measured what it takes to display its power-law tail. Instead, they 
show the mean-squared-displacement of a diffusing micro-sphere. That quantity is 
essentially the velocity auto-correlation function integrated twice, and consequently 
contains the same power-law integrated twice. 

The amplitude that was measured in all these experiments, albeit indirectly with 
photon correlation spectroscopy, is the first-order term in the expansion of P{f) 
above, Eq. (15), in powers of {f /fu)^^'^- This coefficient has two contributions: One 
from the denominator, from Stokes' frequency-dependent friction coefficient, and 
another from the numerator The latter is half-as-large as the former, and with oppo- 
site sign. It stems from the noise term's frequency dependence. 

Instead of measuring a photon correlation function for laser light scattered off 
a suspension of micro-spheres, developments in instrumentation [22] and data anal- 
ysis [21] for optical tweezers have made it possible now to measure directly, with 
accuracy and precision, on a single micro-sphere [24]. Thus it just might be possible 
to observe directly the "color" of the thermal noise, the frequency dependence of 
the non-white power spectrum, in a very challenging single-particle experiment with 
optical tweezers [30]. 
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8 In Situ Calibration of Optical Tweezers by Forced Nano-Scale 
Motion 

There are many ways to calibrate an optical trap. Some ways are better than others 
if accuracy and precision is a concern. In that case, the best way is based on the 
motion's power spectrum [21]. Two aspects must be calibrated: The spring constant 
of the Hookean force exerted by the trap on a trapped micro-sphere (bead), and, 
to this end, the millivolt-to-nanometer calibration factor The latter tells us which 
nanometer-displacement of the bead in the trap corresponds to a measured millivolt- 
change in output potential of a photo diode in the position detection system used 
with the tweezers. A common way to determine this calibration factor requires that 
one knows the radius of the bead, the temperature and dynamic viscosity of the fluid 
surrounding it, and its distance to the nearby surface of the microscope cover slip, if, 
as is usually the case, the experiment is done near this surface. One can then calculate 
the bead's diffusion coefficient in m^/s using Stokes' law (4), Einstein's relation (3), 
and Faxen's formula [31, 32] [21, Sect. XI]. By comparing the result with the same 
quantity measured experimentally in V^/s, the calibration factor is determined. 




frequency (Hz) 



Fig. 2. Power spectrum of 1 .54 /im diameter silica bead held in laser trap with corner frequency 
/c = 538 Hz. The sample moves harmonically with amplitude A — 208 nm and frequency 
/stage = 28 Hz. The power spectrum shown is the average of 48 independent power spectra, 
sampled at frequency /sample ~ 20 kHz. The total sampling time was 79 s, which is six times 
more than we normally would need to calibrate. It was chosen for the sake of illustration, to 
reduce the relative amplitude of the Brownian motion, i.e., the scatter in the spectrum away 
from the spike at 28 Hz. 
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A calibration of the photo diode that is much less dependent on a priori knowl- 
edge, can be achieved by moving the fluid cell with the bead harmonically relatively 
to the laboratory with the optical trap [33]. With a piezo-electric translation stage 
this can be done accurately with an amplitude of order lOOnm and frequency of 
order 30 Hz. In the laboratory system of reference, the fluid flows back and forth 
through the stationary trap with harmonically changing velocity. This gives rise to an 
external force on the trapped bead in (11), a harmonically changing Stokes friction 
force, 

-Fexternal(i) = 7oWstage(<) = 7o27r /staged COs(27r /stage - ^o)) , (17) 

where A and /stage are, respectively, the amplitude and frequency with which the 
stage is driven, and to is its phase. The amplitude A can be chosen so small that the 
forced harmonic motion of the bead in the trap is masked by its Brownian motion, 
when observed in the time domain. Nevertheless, when observed long enough, the 
forced harmonic motion stands out in the power spectrum of the total motion as a 
dominating spike; see Fig. 2. This spike is the dynamic equivalent of the scale bar 
plotted in micrographs: The "power" contained in it is known in m^ because the 
bead's motion in nanometers follows from its equation of motion and the known 
motion of the stage, measured in nanometers. The bead's motion is measured in 
Volts, however, by the photo-detection system, and the Volt-to-meter calibration fac- 
tor depends on the chosen signal amplification, laser intensity, etc. So calibration is 
necessary. It is done by identifying the two values for the power in the spike: The 
measured value in V^ with the known value in m^ [33]^. This method resembles 
an old method of calibration that moves the bead back and forth periodically with 
constant speed, but harmonic motion has a number of technical advantages. One is 
that the precision of power spectral analysis demonstrated in [21] can be maintained, 
while adding the advantage of not having to know the bead's radius, nor its distance 
to a nearby surface, nor the fluid's viscosity and temperature. On the contrary, the 
combination of these parameters that occurs in the expression (3) for the bead's dif- 
fusion coefficient, is determined experimentally from its Brownian motion, so, e.g., 
the bead's radius is measured to the extent the other parameter values are known. 
But also, this calibration method can be used in situ, where an experiment is to be 
done, by confining the bead's forced motion to this environment. This is useful for 
measurements taking place near a surface, in a gel, or inside a cell. 

^ A spike similar to the one shown here in Fig. 2 is seen in [34, Fig. lb]. It was produced with 
a bead embedded in polyacrylamide, hence not moving thermally, and not optically trapped. 
It was used to demonstrate the high sensitivity of the authors' position detection system. 
It was also used for Volts-to-meters calibration of the detection system, and gave 10% 
agreement with the same calibration factor obtained from the power spectrum of Brownian 
motion. The optical properties of polyacrylamide differ from those of water, however, so it 
is an open question how accurate that calibration method can be made. Obviously, it is not 
an in situ calibration method. 
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9 Biological Random Motion 

Robert Brown did not discover Brownian motion, and he, a botanist, got his name as- 
sociated with this physical phenomenon because he in 1827 carefully demonstrated 
what it is not, a manifestation of life, leaving the puzzle of its true origin for oth- 
ers to solve. Brownian motion has been known for as long as the microscope, and 
before the kinetic theory of heat it was natural to assume that "since it moves, it is 
aUve." Brown killed that idea. But after Einstein in 1905 had published his theory 
for Brownian motion, Przibram in 1913 demonstrated that this theory describes also 
the self-propelled random motion of protozoa [35]. By tracking the trajectories x{t) 
of individual protozoa, see Fig. 3, Przibram demonstrated that the net displacement 
x{t) — x{0) averages to zero, while its square satisfies the relationship known for 
Brownian motion, 

(d(t)2) ^ 2nii^Dt , (18) 
where naim is the dimension of the space in which the motion takes place. 



Paramaeciuin, 
1 cm ZGichnuiif; = 0,27 mm Ohjekt. 
Zeit zwischen zwei Punkton 4 Sek. 




Fig. 3. Example of Przibram's motility data, a trajectory of a protozoon, hand-drawn with a 
mechanical tracking device operated in real time with a microscope. A metronome was used 
to mark time on the trajectory every four seconds [35]. 

In Einstein's theory D is the diffusion coefficient, and satisfies his famous re- 
lation (3). Przibram found a value for D which was much larger and much more 
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sensitive to changes in temperature than Einstein's relation states. He used this as 
proof that it was not just Brownian motion that he had observed. 

If Przibram, a biologist, had used a better time resolution by marking out points 
in Fig. 3 more frequently than every four seconds, he might also have gotten ahead 
of the physicists in theoretical developments. But he was drawing by hand, marking 
time to a metronome, so marking points closer to 1 Hz must have been a challenge. 

Fiirth, a physicist at the German university in Prague where Einstein had been a 
professor for 16 months in 1911-12, also studied the motility of protozoa. First he 
repeated Przibram's results, apparently without knowing them [36]. Later he found 
that his data [37] were not described by (18). He consequently considered a random 
walker on a lattice, and gave the walker directional persistence in the form of a 
bias towards stepping in the direction of the step taken previously. By taking the 
continuum limit, he, independently of Ornstein [5, 6], demonstrated that for random 
motion with persistence, (18) is replaced by 

{d{tf) = 2n,,^D{t - P{1 - e-*/^)) , (19) 

where P is called the persistence time, and characterizes the time for which a given 
velocity is "remembered" by the system [37]. 

Ornstein solved (5), since known as the Ornstein-Uhlenbeck (OU) process. Its so- 
lution also gives (19), with P = m/7. The physical meaning of the three terms in the 
OU-process does not apply for cells: Their velocities are measured in micrometers 
per hour, so their inertial mass means absolutely nothing for their motion. Friction 
with the surrounding medium also is irrelevant — the cells are firmly attached to the 
substrate they move on — and it is not thermal forces that accelerate the cells. But as 
a mathematical model the OU-process is the simplest possible of its kind, like the 
harmonic oscillator, the Hydrogen atom, and the Ising model. It also agrees with the 
earliest data. Consequently, the OU-process became the standard model for motility. 
We can write it as 

P^ = -v + {2Dy/''n , (20) 
at 

where each component of 77 is a white noise normalized as in (2) and uncorrected 
with the other components. 

{rjit)) = ; (r/,(t')%(t")) = S,^,Sit' - t") . (21) 

Here 5{t) and 5j^k are, respectively Dirac's and Kronecker's (5-functions, and -qit) is 
assumed uncorrected with v{t') for t > t' . Fiirth's formula (19) is a consequence of 
equations (20) and (21), but follows also from other, similar theories. It was often the 
only aspect of the theory that was compared with experimental data, and with good 
reason, considering the limited quality of data. 

Gail and Boone [38] seem to have been the first to model cell motility with (19). 
They did a time study of fibroblasts from mice by measuring the cells' positions 
every 2.5 hrs. Equation (19) fitted their results fairly well. Since then, cell motility 
data have routinely been fitted with (19). Its agreement with data can be impressive, 
and is usually satisfactory — sometimes helped by the size of experimental error bars 
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and few points at times t that are comparable to P. Data with these properties cannot 
distinguish (19) from other functions that quickly approach 2nii^D{t — P). 

Equation (19) is essentially a double integral of the velocity auto-correlation 
function of the OU-process, where 



= {v{0) • v{t)) 



\t\/P 



(22) 



Experimental results for the velocity auto-correlation function are better suited for 
showing whether the OU-process is a reasonable model for given data. But exper- 
imental results for velocities are calculated as finite differences from time-lapse 
recordings of positions. If the time-lapse is short, precision is low on differences, 
hence on computed velocities. Yet, if the time-lapse is longer, the time resolution 
of the motion is poor. The solution is somewhere in between, compensating for lost 
precision with good statistics. Good statistics was not really achievable till computer- 
aided object- tracking became possible. 




I - 



Fig. 4. Isolated human dermal keratinocytes are motile by nature. If not surrounded by other 
cells, they react as if in a wound: They search for other cells of the same kind with which 
they can connect to form skin. Trajectories are formed from 15min time-lapse photography. 
Trajectories as those shown here in the right panel make up the raw data that are analyzed 
statistically to find a suitable stochastic model for the motility of these keratinocytes. The 
black bar is 0.2 mm long. 
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10 Enter Computers 



We recently wanted to characterize the compatibility of human cells with various sur- 
faces by describing the cells' motility on the various surfaces [39]. Computer-aided 
cell tracking — see Fig. 4 — quickly gave us so much data that we found ourselves in 
a new situation with regards to modelling: We were not limited to showing whether 
or not there is agreement between data and a few consequences of a given model. We 
could investigate the model itself experimentally, measure each term in its defining 
equation, check that their assumed properties are satisfied, and whether together they 
satisfy the equation of motion. 

Furthermore, before we checked the equation of motion, we could check whether 
the data are consistent with various assumptions of symmetry and invariance on 
which the equation of motion is based. We found that the cells behaved in a manner 
consistent with the assumptions that their surroundings are isotropic, homogenous, 
and constant in time. This allowed us to average data over all directions, places, and 
times. This in turn improved the statistics of our investigation of the equation of 
motion [39]. 



11 Tailor-IMade Theory Replaces "One Theory Fits All" 

The theory in (20) states that for a given velocity v the acceleration is a stochastic 
variable with expectation value proportional to v, 

i'iU = -v/P . (23) 

Figure SAB shows that this is also the case for experimental data. 
The theory in (20) states also that 

dv ,dv, dv , , , 1 /n , 

i.e., that this quantity in the OU-process is a white noise with the same speed- 
independent amplitude in both directions: parallel and orthogonal to the velocity. 

Figure 5B shows that experimentally the amplitude of the two components of 
this noise are indeed indistinguishable in the two directions, but the two amplitudes 
are clearly not independent of the speed! Here we see the experimental data reject 
the OU-process as model. The distribution of experimentally measured values of the 
noise also reject the OU-process as model. Figure 5D shows clearly that it is not 
Gaussian, as it is in the OU-process. Apart from that. Fig. 5C shows that the noise is 
uncorrected, like in (21), on the time scale where we have measured it. This result 
radically simplifies the mathematical task of constructing an alternative to the OU- 
process on the basis of experimentally determined properties of these cells' motility 
pattern. 
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Fig. 5. Statistics accumulated from trajectories like those shown in Fig. 4. A: The two com- 
ponents of the acceleration, as functions of speed. Panels Al and A2 show the acceleration 
parallel with, respectively orthogonal to, the velocity. These scatter plots show that the two 
functions contain random parts, like the acceleration in (20). B: Data points with error bars: 
Mean and standard deviation as function of speed for data shown in Panel A. Curves show 
the same quantities, plus/minus one standard deviation, calculated from the theory given in 
(25). C: Correlation functions for scatter shown in Panel A. Panels CI and C2 show the auto- 
correlations of the two components, C3 and C4 show the cross-correlation between the two, 
for both signs of the time difference. The many values shown are almost all indistinguishable 
from zero. This suggests that the scatter in data can be modelled with uncorrelated noise, as 
in (21). This is an experimental result for the theory we seek. The curves shown are not fits to 
the data shown, but results of (25) after it has been fitted to data in Panels B, E, and F. D: His- 
tograms of scatters shown in Panel A, measured relatively to the means shown in Panel B, 
and in units of the standard deviations shown in Panel B. The curves shown are not fits to the 
histograms shown, but results of (25) after it has been fitted to data in Panels B, E, and F. 
E: Velocity auto-correlation function, calculated from trajectories like those shown in Fig. 4. 
It is not a simple exponential as in (22). But a sum of two exponentials fit data perfectly. So 
we assume that the theory we seek has a velocity auto-coiTelation function that is a sum of two 
exponentials. The curves through the data points are that correlation function, plus/minus one 
standard deviation, computed with the theory in (25), after it has been fitted to the data shown 
here, and simultaneously to the data in Panels B and F. F: Histograms of speeds and (speed)^ 
read off trajectories like those in Fig. 4. The curves shown are the same speed distributions 
calculated from the theory in (25), after it has been fitted to the data. 



The velocity auto-correlation function of the OU -process is a simple exponential, 
(22). Figure 5E shows the experimentally measured velocity auto-correlation func- 
tion. It is fitted perfectly by the sum of two exponentials, so again the experimental 
data reject the OU-process as model. 

The data shown in Fig. 5 are so rich in information that with a few assumptions 
favored by Occam's Razor one can deduce from the data which theory it takes to 
describe the data, and this theory is unambiguously defined by the data [39]. Results 
from this theory are shown as the fully drawn curves passing through the data points 
in Fig. 5. It is given by the stochastic integro-differential equation 

dii 

— it) - ^iSvit) (25) 
W f dt'e-"''''-*'^v{t')+a{v{t))r]{t) , 

where 

a{v) = (To + f^iv ■ (26) 

The integral over past velocities in (25) is called a memory-kernel by mathemati- 
cians. It shows that these cells have memory. This is no surprise: The polarity of 
the cytoskeleton of a moving cell is a manifest memory of direction, and while its 
instantaneous velocity depends on the activity of transient pseudopodia, the fact that 



Brownian Motion after Einstein 17 



pseudopodia are active depends on states of the cell that last longer than the individ- 
ual pseudopod, one would expect. 

Note the similarity between (25) and the hydrodynamically correct theory for 
Brownian motion, (10,11). Though both are more complex than Einstein's theory, 
applied by Przibram, and the OU-process, applied by Fiirth, they still have much in 
common. This is so because they both are linear and both respect causality and the 
same space-time symmetries. 

Linearity is simplicity, so wherever in modelling it is sufficient, one avoids going 
beyond it. This is why (10,1 1) and (25) are both linear. 

The Principle of Causality states that the future does not affect the present, in- 
cluding present rates of change of state variables. Only the past can do this. This prin- 
ciple is respected throughout physics, and we have of course built it into our motility 
models as well. This is why the rate of change of the velocity given in (10,11), re- 
spectively (25), depends only on past and present velocities. The integral kernels 
occurring in both equations are memory kernels in order to respect this principle. 

In a homogenous, isotropic environment that is constant in time, there is no ab- 
solute position, direction, nor time. A theory for a dynamical system in such an en- 
vironment consequently cannot depend on the position variable x, nor can it depend 
explicitly on the time variable t, nor on explicit directions in space. The theory must 
be translation invariant in space, time, and with respect to direction. The last invari- 
ance is called covariance under rotations, because a theory for a vector variable like 
the velocity is not invariant under rotations of the coordinate system, it is covariant, 
i.e., transforms like the vector it describes. Because these space-time symmetries are 
shared by hydrodynamics and our cells, neither (10,1 1) nor (25) depends on x, nor 
explicitly on t, and both models transform like a vector under rotations. 

We conclude that with the rich data that one now can record and process, one 
should not be satisfied with the simplest possible model for persistent random mo- 
tion, the OU-process. "One size fits all" is no longer true, if it ever was. Motility 
models can be made to measure. Here we have only presented the first phenomeno- 
logical steps of that process: How to plot and read motility data in a manner that 
reveals mathematical properties of the theory sought. That done, it is another task to 
construct a model with the properties demanded. If that can be done, it is yet another 
task to decide whether the theory is unique or not. Two examples of such theories 
and their derivation are given in [39]. 
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